Vacuum fluctuations and generalized boundary conditions 
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We present a study of the static and dynamical Casimir effects for a quantum field theory satisfying 
generalized Robin boundary condition, of a kind that arises naturally within the context of quantum 
circuits. Since those conditions may also be relevant to measurements of the dynamical Casimir 
effect, we evaluate their role in the concrete example of a real scalar field in 1 + 1 dimensions, a 
system which has a well-known mechanical analogue involving a loaded string. 
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I. INTRODUCTION 

Over the course of the last 15 years, there has been a 
renewed interest in the Casimir effect. This was partly 
due to a second generation of experiments, started in 
1997, which triggered a sustained flow of both theoretical 
and experimental works. As an outcome of those works, 
our knowledge about the dependence of Casimir forces 
on the geometry and material properties of the objects 
involved has remarkably improved jl| . 

The dynamical counterpart of the Casimir effect, also 
known as 'motion induced radiation', has also been the 
subject of intense research, manifested in the consequent 
profusion of works 0] ■ Although the direct measurement 
of radiation generated by moving mirrors is a daunt- 
ing experimental challenge, photon creation induced by 
time dependent external boundary conditions has, in- 
deed, been observed experimentally, albeit in a different 
context, namely, superconducting circuits Q. There are 
also ongoing experiments aimed at measuring the photon 
creation induced by the time dependent conductivity of a 
semiconductor slab enclosed by an electromagnetic cav- 
ity [if , and proposals based on the use of high frequency 
resonators to produce the photons, and ultracold atoms 
to detect the created photons via superradiance Q . 

The superconducting circuit experiment mentioned 
above Q, consists of a coplanar waveguide terminated 
by a SQUID, upon which a time dependent magnetic flux 
is applied. This system can be described by a quantum 
scalar field <f>(x, t) (the magnetic flux at the different po- 
sitions of the transmission line) satisfying on the SQUID 
(located at x = a) the boundary condition 



L dx 



o2 ± 

{a,t) + C-^{a,t) + E{t)4>{a,t)=Q, (1) 



where C and L are constants, and E(t) is proportional 
to the (time-dependent) external magnetic flux. Equa- 
tion JT]) can be interpreted as a sort of generalized Robin 
boundary condition, because of the presence of the term 
with second-order time derivatives. The theoretical anal- 
ysis of that experiment was done assuming that the term 



proportional to the second derivative of the field is neg- 
ligible @. Whence the boundary condition became a 
standard Robin boundary condition with a time depen- 
dent, externally driven parameter. The boundary condi- 
tion that results after implementing this approximation, 
has been described in terms of an effective length 0, Q . 

From a theoretical point of view, the particle creation 
rate for the case of a quantum field satisfying (time in- 
dependent) Robin boundary conditions on a moving mir- 
ror has been investigated in Ref. @- The complemen- 
tary situation, namely, time dependent Robin boundary 
conditions on a static mirror have been considered in 
Refs. [1, It is the aim of this work to consider the 
static and dynamical Casimir effects for quantum fields 
which are subject to generalized boundary conditions of 
the type defined in Eq.Q. Our interest in this prob- 
lem is twofold. On the one hand, the presence of the 
second derivative of the field in the boundary condition 
modifies the usual Sturm-Liouville problem, what man- 
ifests itself in the existence of an eigenvalue-dependent 
boundary condition. This problem, well-known to math- 
ematicians, has also been considered in different areas 
of physics. It has not, however, to our knowledge, been 
dealt with in the Casimir effect context. 

Regarding the experimental relevance of the inclusion 
of this kind of term, we note that, although it may be ne- 
glected in the particular experimental setup considered in 
Ref. Q , this contribution to the boundary condition may 
indeed affect the spectral distribution of the created par- 
ticles [sl, something to which future experiments might 
be sensitive. 

This paper is organized as follows. In Section (TTJ we de- 
scribe the classical aspects of the model, including a me- 
chanical analogue, a discussion of the eigenvalue problem 
and the necessity of modifying the inner product between 
cigenfunctions due to the presence of the second order 
derivative in the boundary condition. In Section IIIII we 
compute the static vacuum energy, Section IWl deals with 
the dynamical Casimir effect, and Section |V] contains our 
conclusions. 
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II. CLASSICAL MODEL: THE LOADED 
STRING ANALOGY 

A. A mechanical analogue 

Many different boundary conditions a scalar field in 
f + f dimensions can be subject to, can be realized in 
classical mechanics analogues, based on vibrating strings. 
Indeed, Dirichlct boundary conditions correspond to a 
string with fixed endpoints, while Neumann conditions 
can be implemented by attaching, to the corresponding 
cndpoint, a massless ring which can slide freely and fric- 
tionless along a vertical rod. Robin boundary conditions 
result, in turn, when the ring is also coupled to a vertical 
spring @. Finally, the generalized boundary condition 
of Eq.((T]) can be generated by letting the ring be mas- 
sive rather than massless (in all cases we will consider 
the linear regime, i.e. small amplitude oscillations of the 
string). 

To see this, let us assume T to be the string tension, 
\x its mass density, and that its configuration may be de- 
scribed by a single function y{x,t), measuring its vertical 
departure from the equilibrium configuration. Denoting 
by m the mass of the ring and by k the spring constant, 
the position of the ring, y(a,t) satisfies Newton's equa- 
tion: 



my(a, t) = -ny(a, t) - Ty'{a, t) 



(2) 



which has the same form as the generalized boundary 
condition of Eq.([T]). The vibrating string problem with 
this kind of condition on one of its endpoints constitutes 
a well known problem in classical vibrations. 

We stress that the boundary condition for the defor- 
mation of the string is in fact the dynamical equation for 
the position of the ring q(t) = y(a, i), and this is the ori- 
gin of the presence of second order time derivatives in the 
generalized boundary condition. A proper treatment of 
the system should regard q(t) and y(x,t) as qualitatively 
different degrees of freedom, in the sense that q has a 
discrete, finite mass, while y is endowed with a continu- 
ous mass density. An enlightening way to treat this kind 
of problem can be seen, for example, in (Toj . where the 
quantization of a non relativistic string with an arbitrary 
mass distribution has been considered. There, a system 
with a single continuous mass distribution which approx- 
imates the mixed continuous and discrete case has been 
considered. In this way, all the steps of the Lagrangian 
and Hamiltonian formalisms, and even canonical quan- 
tization, are well defined. The desired mass distribu- 
tion is approached at the end of the process, as a special 
limit, after all the stumbling blocks in the procedure are 
avoided. 

We will follow here an alternative procedure, which, 
as we have explicitly checked, yields the same results. 
Since the difficulties in the problem at hand come from 
the fact that the discrete mass is precisely at one of the 
endpoints of the system, we avoid the coincidence of those 
two singularities by temporarily splitting them. Indeed, 



we shall first assume that the mass m is located at an 
arbitrary position xq, < xq < a, and impose Neumann 
boundary conditions at x — and x = a. The generalized 
Robin boundary condition at x = a is then recovered by 
taking the 'coincidence limit' xo —> a~ . 
The classical Lagrangian then reads: 



1 f a 

L = — I dx 

2 Jo 



-Ky S(x - x ) + m ( — ) 6(x - xq) 



T 



dy 



dy 
dx 

2 



, (3) 



with Neumann conditions implicitly assumed at x = 
and x = a. 

From the classical equation of motion one can easily 
check that the presence of a localized mass on the string 
induces a discontinuity in the spatial derivative of y, 

T ~ = m v( x o,t) + Ky(x ,t) . (4) 

Note that the string can interchange energy with the 
mass, and the conserved total energy of this system reads 



2 



dx 



l fltJ +T (dx 



-Ky 2 (x ,t) 



1 (dy 

2 \dt 



(5) 



that is, the sum of the mechanical energies associated to 
the string and the ring. 



B. A one dimensional cavity with localized 
conductivity and permittivity 

Let us now focus on the analogous case of a scalar field 
in 1 + 1 dimensions, as a toy model for the electromagnetic 
field in 3 + 1 dimensions. We assume the Lagrangian to 
be given by the expression: 



2 



dx 



-<A ( d(j> 



with 



and 



-Vix,^ 2 



e(x, t) = 1 + \(t)S(x - Xq) 



V(x,t) = v(t)S(x - Xq) 



(6) 



(7) 



(8) 



As in the mechanical model, we shall regard this as a sim- 
ple model to describe a cavity in which the permittivity 
and conductivity are concentrated at the point x = x , 
which, tending to a from the left, reproduces the gen- 
eralized Robin condition at x — a. The particular case 
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X(t) = and v(t) — v (l + f(t)) has been considered in 
Ref. 11[ , as a simple model of the experimental setup of 
Ref. 4( (the generalization to the electromagnetic case has 
been analyzed in Ref. (l2j ) . In Section Hill we will com- 
pute the Casimir interaction energy between two slabs 
described by constant values of v = vq and x = Xo- In 
Section IV we will compute the photon creation associ- 
ated to a time dependent v(t) and x — Xo- 



C. Eigenfunctions and inner product 

Let us consider the model given in Eq.([B]). with Neu- 
mann boundary conditions at x = 0, a. We will compute 
the eigenfunctions and eigenvalues for the particular case 
v(t) = Vo and e(t) = Xo- The eigenmodes can be written 
as 

^ k (x,t) = N k e~ tkt [cos(fcx) cos(fc(x - a))0(x o - x) 



+ cos(/ccco) cos(k(x — a))9(x — xq)] 



-ikt 



M*), (9) 



where Nk is a normalization constant. They are continu- 
ous at x = xq and satisfy Neumann boundary conditions 
at x = 0, a. The discontinuity equation at x = xq implies 

k sin(fca) = (vq — Xok 2 ) cos(fcxo) cos(fc(xo — a)) , (10) 

which is the equation that defines the eigenfrequencies. 
In the particular case xq — » a, the transcendental equa- 
tion that defines the eigenfrequencies simplifies to 



fctan(fca) = (vq — Xok 2 ) . 



(11) 



It is straightforward to show that, unless xo = 0, eigen- 
functions corresponding to different eigenvalues are not 
orthogonal with the usual inner product. However, defin- 
ing a generalized inner product (see the Appendix and 
Refs. [ll EDI) 



(ipk ,i>k>) = dx ip k (x)tpk' (x) + xo^Pk (a) ipk' (a) 
Jo 



dx e(x)ip k (x)ipk' (x) 



(12) 



one can check the orthogonality (tpk, i'k 1 ) = for k ^ k'. 
With appropriate normalization, the eigenfunctions can 
be chosen to be orthonormal, as we shall assume it has 
been done in what follows. 

This phenomenon is a general feature of the hybrid 
continuous plus discrete systems, where one can show 
that the eigenfunctions of the Hamiltonian are orthogonal 
for a scalar product defined in terms of a kernel fiol ]. 
which defines a Sturm-Liouville problem. 

The equations above are valid even in the time- 
dependent case wo - > v (t)i m which the eigenvalues be- 
come paramctrically dependent on time, as well as the 
eigenfunctions ipk- We will analyze the time-dependent 
situation in Section IV. 



Writing the field 
tial eigenfunctions 



as a linear combination of the spa- 



{x,t) = ^2Q k (t)i>k(x).. 



(13) 



and inserting this expression into the classical Lagrangian 
one can check that it reduces to a set of uncoupled har- 
monic oscillators Qk, with frequency k: 



L 



2 ^ 

h 



(Ql-k 2 Ql) 



(14) 



Details of this calculation are presented in the Appendix. 
Note that the additional term in the inner product 
Eq. ([T2"]) is crucial to cancel the kinetic term concentrated 
at xq. Note also that these results can be straightfor- 
wardly generalized to cases where several slabs are lo- 
cated at different positions. 



III. STATIC CASIMIR EFFECT 

Once the system has been reduced to a set of uncou- 
pled harmonic oscillators, the calculation of the vacuum 
energy can be performed by direct mode-summation. We 
will explicitly compute the static Casimir effect for two 
different physical situations: the interaction between two 
thin slabs, each one described by its conductivity and 
permittivity (in free space), and a system that satisfies 
generalized boundary conditions. 



A. Interaction vacuum energy for two slabs 

Let us consider two slabs, located at x = ±a/2, and 
described by permittivities Xo" an d conductivites v ± , re- 
spectively. In order to have a discrete set of eigenfre- 
quencies, we enclose the system in a box of size 2L, and 
impose Neumann boundary conditions at x = ±L. As we 
will take the limit L 3> a at the end of the calculation, 
the result will be independent of L and of the boundary 
condition imposed at x = ±i. 

The eigenfunctions that satisfy Neumann boundary 
conditions at x — ±X can be written as 

!Ai cos[cj(.t + L)} for — L < x < — f 
A 2 cos(wx) + A 3 sin(wx) for — | < x < | 
Ai cos[cj(.t - L)] for f < x < L . 

(15) 

The function f u {x) must be continuous at x = ±a/2 and 
the spatial derivatives must satisfy: 

disc[d^U ±Q/2 = [v± - Xo" 2 \U±\) ■ (16) 

The eigenfrequencies are the solutions of det M = 0, 
where M(a,L,u>) is the 4x4 matrix associated to the 
linear system of equations for the coefficients Ai derived 
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by inserting Eq. (fT5|) in Eq. (fT6l) . After some straightfor- 
ward calculations one can show that 

det[M(a,L,u)} = A+A~(sin(2w(a - L)) + 2sin(aw)) + 

2(A+ + A") cos(aw) + (A+A~ - 4) sin(2Lw) + 

2(A+ + A")cos(2iw) , (17) 

where A ± = ±(v± -Xo^ 2 )- 

We will compute the Casimir energy as the difference 
between the zero point energy of the slabs separated by 
a distance a, and that corresponding to a distance / 3> a. 
Using the argument theorem, we see that: 



E c (a) 



2 ^ 



1 

47TI 



= — - — r f dz logdet 



M(a, L, z) 
M(l,L,z) 



(18) 



where oj n and w n are the eigcnfrequencies associated to 
the distances a and I, respectively. The integration path 
must include the real positive axis. Following standard 
steps, and taking the limit L — » oo, we arrive to an inte- 
gral in the imaginary- frequency axis z = i£: 



E c = 



/ d£ log [1 - e- 2 « a C(0] , (19) 

J<3 



The eigcnfrequencies are defined implicitly by Eq. (TTT]) . 
that we rewrite as G(a,u>) = 0, with 

G(a, uj) — u> sin(a;a) — (vq — Xow 2 ) cos(wa) . (21) 

As in the previous subsection, we enclose the system in 
a large box of size 2L, and impose Neumann boundary 
conditions at x = ±L. The eigcnfrequencies are thus 
determined by the equation G(a, L, u) = 0, with 




where 



c(0 



xoe) 



(20) 



Not surprisingly, for the particular case \o =0, this 
result coincides with the usual Casimir energy for the so 
called J-potentials. Moreover, in the Vq — > oo limit, one 
has C(^) — > 1, and the result reproduces the usual one for 
Dirichlet boundary conditions. Note that, as dC/dxt > 
0, the presence of the second order time derivative in the 
boundary conditions enhances the interaction between 
slabs. 

It is interesting to remark that the final result for the 
Casimir energy is tantamount to the one corresponding 
to a (5-potential with a frequency-dependent coefficient 
v(u) = vo — Xo 1 ^ 2 - Therefore, this result could have been 
derived using Lifshitz formula with the particular reflec- 
tion coefficients that describe the slabs. The above cal- 
culation is an alternative and equivalent way to compute 
the vacuum energy, that shows that the Casimir energy 
is just the sum over the eigenfrequencies defined by the 
w-dependent boundary conditions. 



B. Interaction vacuum energy for generalized 
boundary conditions 

We now consider the second case, namely, a system 
that satisfies Neumann boundary conditions at x = 
and generalized Robin boundary conditions at x = a. 



FIG. 1. Dimensionless total Casimir energy Ec/vo as a func- 
tion of the distance voa, for different values of voxo- While the 
force is always repulsive for large values of voXo, for voXo — 0.1 
(solid line) it starts being repulsive, then becomes attractive, 
and is finally repulsive again for larger values of the distance. 
Dashed line corresponds to VoXo = 0.5, and the dotted line 
corresponds to i>oXo ~ 5 



G(a, L, lo) = G(a, ui) sin[w(L - a/2)f 



(22) 



Then we compute the Casimir energy as the difference 
between the vacuum energy associated to the length a, 
and that associated to I with I ;§> a. Using again the 
argument theorem 



E c {a) 



1 

47TI 



dz log 



G(a, L, z) 
G(l,L,z) 



the final result can be written as 
1 



E, 



c 



2- 



d£ log(l - e- 2 ^ a B(0) , 



(23) 



(24) 



o 



where the function B contains the information about the 
generalized boundary condition 



B(0 



£ + v + xoe 



(25) 



As the sign of B depend on the values of vq, xo and £, 
the force can be attractive or repulsive depending on the 
value of a. 
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In oder to analyze the behavior of the energy with the 
different parameters, it is useful to note that 

E c 

= h(v a, voxo) , (26) 

«o 

where h is a dimensionless function. This can be eas- 
ily checked changing variables £ — > £,/vq in Eq . (|24| . 
Therefore, the qualitative dependence of the energy with 
the distance only depends on the dimensionless quantity 
«oXo- For instance, when vqXo > 1/4, B is negative 
for all values of £, and therefore the force is repulsive. 
In particular, in the limit vq — > oo, the boundary con- 
dition at x = a becomes Dirichlet boundary condition, 
and one has B — > — 1. The Casimir energy becomes 
the standard result for a scalar field satisfying Neumann 
boundary conditions at x = and Dirichlet boundary 
conditions at x = a, which corresponds to a repulsive 
force. When VoXo < 1/4, we cannot predict the sign of 
the force analytically. 

In Fig. 1 we present some numerical evaluations of 
Eq. ([2^|) . We refer all the ingredients in the expression 
for the energy to the dimensionful quantity, vo. We plot 
the Casimir dimensionless energy Eq/vq as a function 
of the distance a (in units of vq) for different values of 
voXo- We see that for small distances, the force changes 
sign for the smaller value of voXo = 0.1. In this regime, 
the force cross from repulsive to attractive and back sign 
to repulsive again, as the distance increases. For other 
values of vqXo > 1/4 the force is, as expected, always 
repulsive. 

IV. DYNAMICAL CASIMIR EFFECT 

Let us now consider the case of a time dependent v(t) = 
vq(1 + f(t)) and constant x(t) = Xo- In order to analyze 
this problem one can proceed as usual @, [ll|, 13 s showing 
that, at the classical level, the system can be described by 
a set of coupled harmonic oscillators with time dependent 
frequencies and couplings. To this end, we introduce an 
'instantaneous basis' il>k(x,t) through the equations 

Vi' + fe 2 (t)Vfc = o 
V4(o) = o 

1^(a) = 

Vi(4) - Vi(^) = (-Xofc 2 + v(W k . (27) 

Note that, due to the time dependent boundary condi- 
tion at Xo, the eigenvalues are time dependent, and this 
induces a parametric time dependence in the basis func- 
tions. The limit xo —> a~ describes a cavity with a 
SQUID at one end, but one could consider more gen- 
eral situations like a cavity ended by two SQUIDs, or 
even a set of SQUIDs located at different positions in a 
waveguide. 

We now expand the field using the instantaneous basis 
<(>(x,t) =^2Q k (t)Mx,t). (28) 

k 



The classical Lagrangian given in Eq.© can be writ- 
ten in terms of the variables Qk, as in the static case. 
The frequencies of the classical oscillators become time- 
dependent, and the Lagrangian contains additional terms 
proportional to derivatives of the basis functions, which 
in turn are proportional to the derivatives of the eigen- 
values k{t). 

As shown in the Appendix, the classical Lagrangian 
reads 

L = \ £Wfc - fc2 (t)Ql) + E A kj(t)QkQj 

k kj 

-\Y. S ^)QkQj, (29) 

kj 

where the time-dependent matrices Akj and Skj can be 
chosen to be antisymmetric and symmetric, respectively. 
The explicit expressions for them are derived in the Ap- 
pendix. Note that, if the time dependence of the eigen- 
values is proportional to a dimensionless parameter 5, 
A kj = 0(6) and S kj = 0(5 2 ) (sec Eq.^TBJ). 

Let us now assume that the externally driven prop- 
erty v(t) has a harmonic time dependence with frequency 

0. When the external frequency is tuned with an eigen- 
frequency of the static cavity f2 = 2fco, with fco one of 
the solutions to Eq. (fT0|) . we expect the number of cre- 
ated photons to be enhanced by parametric resonance. 
Moreover, since in general the spectrum of the cavity 
is not regularly-spaced, it is reasonable to neglect cou- 
plings between modes [l6j. The dynamics of the system 
is essentially given by that of the mode associated to fco, 

1. e. a single harmonic oscillator with time dependent fre- 
quency: 

L^l[Ql-k 2 (t)Qt Q ], (30) 

where we have neglected terms of order 0(5 2 ). Here k (t) 
is a solution to 

fcsin(fca) = (v(t) — Xok 2 ) cos(fccco) cos(fc(xo — a)) . (31) 

For the sake of simplicity, we will solve this equation 
for xo — >• a~ . We denote by fco a solution to Eq. (|3"Tj) in the 
static case f(t) = 0. Assuming that ko(t) = fco (1 +?//(£)), 
with 7]f(t) -C 1, one can show that 



V + kl{a + xo) + a{v - Xo^o) 2 

When f(t) = Asin[2kot], the number of photons with 
frequency fco will grow exponentially. The calculation is, 
by now, standard 0, [TJ [T3, [l6| , and will not be repro- 
duced here. The result is 

A^ o (t) = cxp[fc ^] = exp[Ai] . (33) 

Note above equations are valid for small values of the 
aforementioned parameter 8 = At] <C 1. 
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It is interesting to analyze two opposite limiting cases. 
When vqCl ^> 1, the lowest frequency solution is k ~ 
7r/(2a), and therefore r\ ~ l/(voa). The particle creation 
rate, in this case, is A = irA/(2voa 2 ). 

On the other hand, when vqo, <C 1, the first solution to 
the transcendental Eq. (fTT|) is 



vpa 



(34) 



In this case we have 77 « 1/2 (assuming that Xo/o *C 1). 
Therefore, the corresponding rate is A = A/2^/vo/ a. 

The last case may be of some interest for the experi- 
mental observation of the dynamical Casimir effect, since 
the Robin or generalized Robin boundary conditions may 
be adjusted in such a way to reduce the value of the low- 
est eigenfrequency of the unperturbed cavity. This point 
deserves further analysis. 

When considering the parametric resonance situation, 
it is crucial to tune the external frequency with one of the 
eigenfrequencies of the system. The term proportional 
to xo m the boundary condition, however small, does 
introduce significant modifications to the static eigenfre- 
quencies of the cavity. This effect may be relevant, for 
instance, for an experiment with a coplanar waveguide of 
finite size. 

As an illustration of the last point, let us assume that 
xq — > a~ and that the external frequency is twice the 
lowest eigenfrequency of the cavity (that is, twice the 
first solution to Eg . (Til"]) ) . In a realistic situation [l5| . 
VqCl 3> 1, and an approximate solution is fc — 7r/(2a). 
Expanding the transcendental equation around this so- 
lution we obtain 



2a 



(1 



I 

v a 



An 



2„2 



) 



(35) 



Assuming [I5| voa — 20 and Xo/a = 0.05 the correc- 
tion to the lowest eigenfrequency due to xo is 3 x 10 ~ 4 . 
This small correction should be taken into account in 
order to achieve parametric resonance if the amplitude 
of the time-dependent external conditions is sufficiently 
small (H. 



V. CONCLUSIONS 

The original aim of this work has been to analyze the 
static and dynamical Casimir effects when the field sat- 
isfies generalized boundary conditions involving second 
order time derivatives. These conditions were in turn 
motivated by the effective boundary condition satisfied 
by the magnetic flux in a waveguide terminated by a 
SQUID, a setup that has been recently employed as a 
device to measure the creation of photons from the vac- 
uum in the presence of time-dependent external fields Q . 

We have shown that this problem does have a sim- 
ple and well-known classical analogue: a loaded string. 
From this point of view, the presence of the second-order 



time derivative is not surprising, since the boundary con- 
dition is in this case nothing but the dynamical equa- 
tion for a massive ring attached to the end of the string. 
From a mathematical point of view, we have eigenvalue- 
dependent boundary conditions, and therefore the eigen- 
functions associated to different eigenvalues are not or- 
thogonal under the usual inner product. A generalization 
of the inner product makes them orthogonal. When ex- 
panding the deformation of the string in terms of spatial 
eigenfunctions, the additional term in the inner product 
exactly cancels out the kinetic term associated to the 
ring, and one ends up with a set of decoupled harmonic 
oscillators. 

This mechanical analogy lead us to consider a more 
general situation, in which the ring is not attached to the 
end of the string; rather, it is located at an arbitrary dis- 
tance from the endpoint. From a field-theory perspective, 
this can be considered as a toy model for an electromag- 
netic cavity in which one inserts a thin slab characterized 
by its conductivity and permittivity. From the quantum 
circuits point of view, this corresponds to a situation in 
which a SQUID is inserted in a one-dimensional waveg- 
uide. 

Therefore, we computed the static Casimir interaction 
energy between two slabs, generalizing previous results 
for ^-potentials. We also computed the static vacuum 
energy for the particular case in which the slab is near a 
border of the cavity. In summary, we have shown that the 
Casimir energy can be computed as the sum over modes 
satisfying the generalized Robin boundary condition. 

Finally, and coming back to the original motivation, 
we considered some particular aspects of the dynamical 
Casimir effect for the generalized boundary conditions. 
On the one hand, we computed the particle creation rate 
assuming parametric resonance for the case of a finite 
waveguide ended by a SQUID. On the other hand, we 
discussed the influence of the second order time derivative 
on the tuning of the external pumping with one of the 
eigenfrequencies of the cavity. We have seen that the 
term proportional to </> in the boundary condition may 
indeed be relevant for this tuning. 

The analysis of loaded strings with masses distributed 
periodically along it, or analogue acoustic and elastic sys- 
tems, induces the presence of (approximate) band gaps 
in the spectrum [17|. It would be interesting to ana- 
lyze theoretically electromagnetic analogues of these con- 
figurations, like a waveguide with several SQUIDs, dis- 
tributed periodically on it, or a microwave cavity with 
several slabs inserted accordingly. 



Appendix A 

In this Appendix we present some details of the calcu- 
lations of the classical Lagrangian in both the static and 
dynamical cases. 

To begin with, let us show that it is necessary to modify 
the usual inner product in order to have an orthonormal 
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basis. The eigcnfunctions satisfy 

i>' k ' + k 2 i, k = 
with the boundary conditions 



Therefore 



(Al) 



^(0) = o 
^(a) = 

^(4) - 4i(*o) = (-Xofc 2 + «o)Vfc ■ (A2) 

For the sake of definiteness, we choose Neumann bound- 
ary conditions at x = 0, a. The results below can be 
generalized to the case of Robin boundary conditions. 
We compute 



hj= / m j -i/) i il/J)dx = i3-n / Mjdx, (A3) 
Jo Jo 

where in the last equality we used the eigenvalue equa- 
tion. On the other hand we have 



One should be careful with the evaluation of 



(A4) 



(A5) 



because of the discontinuity at xq. Eq. (|A4[) can be ob- 
tained writing J Q = J Q + and using the boundary 
condition at xq: 



= (Xof ~ vo)il) k (xo)%l)j(xQ) (A6) 



(A7) 



Subtracting Eqs. (|A3j) and (|A4[) we get 

0=(J 2 -^ 2 ) 
and therefore the inner product defined as 



o 



{<t>i,4>j) = xoi>i(xo)ipj(xo) + / dxijjiipj 

Jo 



dx e(x)ipi^j 



(A8) 



vanishes when i ^ j. 

To compute the static classical Lagrangian we write 



4>(x,t) = J2Qk(t)Mx)- 



(A9) 



pa pa 

/ dxcj) 2 = y^Q k Qj I dxip k i(jj 
Jo kj Jo 



^2QkQj{Skj - Xotpk{xo)ipj{x )) 



kj 



(A10) 



A similar calculation can be done for the spatial deriva- 
tives 



f dx<f>' 2 = V Q k Q f dxiP'^ 
Jo kj Jo 

= ^2Q k Q 3 [ a [(ipkW -TpkTp'-)} 



kj 



kj 



J^QkQ, / [(Vfc^O' + iVfc^OKAii) 



Inserting Eq. (|A6J) in (|Allj) . and using again the orthog- 
onality 

f ' dx<p' 2 = yVQ 2 -^ 2 (*<),*)■ (A12) 
Jo k 

Using Eqs. ()Al0p and (|A12[I . the classical Lagrangian 
Eq.©, can be written in the static case as 



(A13) 



We now consider the time dependent situation v — > 
v(t). Using the basis functions introduced in Section IV, 
from Eq.tpg]) we have 

4(x,t) = Y^(Qk(t)ij>k(x,t) + Qk(t)il>k(x,t)), (A14) 



and therefore 

1 r 

2 



where 



f a 1 

/ dx e(x)<j> 2 = - V Ql + V QkQjA kj 

JO * ,, UA 



kj 
kj 



77 QkQjSk 



(A15) 



A k j = / dx e(x) tpkipj 
Jo 



Skj 



dx e(x) i>kipj 



(A16) 



Note that, due to the orthogonality of the eigenfunctions, 
the matrix A k j is antisymmetric. The matrix S k j is ob- 
viously symmetric. 
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